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Abstract. Random sequential adsorption (RSA) models have been studied 
[1] due to their relevance to deposition processes on surfaces. The depositing 
particles arc represented by hard-core extended objects; they are not allowed 
to overlap. Numerical Monte Carlo studies and analytical considerations are 
reported for ID and 2D models of multilayer adsorption processes. Depo- 
sition without screening is investigated, in certain models the density may 
actually increase away from the substrate. Analytical studies of the late stage 
coverage behavior show the crossover from exponential time dependence for 
the lattice case to the power law behavior in the continuum deposition. 2D 
lattice and continuum simulations rule out some "exact" conjectures for the 
jamming coverage. For the deposition of dimers on a ID lattice with dif- 
fusional relaxation we find that the limiting coverage (100%) is approached 
according to the ~ \j \ft power-law preceded, for fast diffusion, by the mean- 
field crossover regime with the intermediate ~ 1/t behavior. In case of fc-mer 
deposition (k > 3) with diffusion the void fraction decreases according to the 
power-law t~ 1 ^ k ~ 1 \ In the case of RSA of lattice hard squares in 2D with 
diffusional relaxation the approach to the full coverage is ~ t~ x /' 2 . 



1. Kinetics of multilayer adsorption 

In this section we study [1] the deposition of fc-mers of length k on the linear 
periodic ID lattice of spacing 1, and the deposition of square-shaped (k x k)- 
mers on the periodic square lattice of unit spacing in 2D. The deposition site 
is chosen at random. The time scale, t, is fixed by having exactly L deposition 
attempts per unit time. We studied systems of sizes up to L = 10 5 in ID and 
L x L — 1000 2 in 2D. Various Monte Carlo runs (for k = 2,3, 4, 5, 10 in ID, 
and k = 2,4 in 2D) went up to t = 150. Our numerical values are consistent 
with exact results for layer n = 1 in ID [2]. 

We first report on a model with decreasing layer coverage. Only if all the 
lattice segments in the selected landing site are already covered by exactly 
(n — 1) layers, the arriving fc-mer is deposited, increasing the coverage to n 
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(n > 1). For small t, the coverage (fraction of the total volume covered by 
depositing particles) increases according to n (t) oc t n , as expected from the 
mean-field theory. Our results suggest that for lattice models the fraction 
of the occupied area in the nth layer, 6 n (t), approaches the saturation value 
exponentially 9 n (t) w 8 n (oo) + B n exp(— t/r n ). However, the jammed state 
in the higher layers in the deposition without overhangs contains more gaps 
the larger is the n value. The growth in the higher layers proceeds more 
and more via uncorrelated "towers". We find that the jamming coverages 
vary according to a power law, with no length scale, reminiscent of critical 
phenomena, 

^(ooJ-Moo) 5 *^ (!) 

Within the limits of the numerical accuracy the values of the exponent are 
universal for k > 2. Based on numerical data analysis we find the estimates 
0(1D) = 0.58 ± 0.08 and 0(2D) = 0.48 ± 0.06. These values are most likely 
exactly 1/2 as suggested by analytical random- walk arguments [3]. 

Now we report on a model with decreasing layer coverage. For layers n > 1 
the deposition is successful only if no gaps of size k or larger are covered. 
Thus, the deposition is always allowed if all the "supporting" k sites in the 
(n— l)st layer are filled or have only small internal gaps. The coverage at short 
times decreases with layer number. However, for the particular deposition 
rule considered here the coverage in layer n eventually exceeds that in layer 
(n — 1) at larger times. This unexpected behavior was found numerically 
for all layers n < 55 and for all k (k — 2,3,4,5,10) studied. We found 
clear evidence of the power-law behavior (eq.(l)), with A < 0. We found 
the power cf) to be universal for all k studied, cf) — 0.3 ± 0.15. When large 
enough covered (by fc-mers or gaps of sizes up to k — 1) regions have formed 
in layer (n — 1), then the deposition with overhangs beyond those regions 
will be delayed. Thus, there will be some preference for higher density in 
layer n especially near the ends of the regions occupied in layer (n — 1). 
To test the above suggestion, wc considered the following monolayer dimer- 
deposition model. We select randomly pL/2 dimers and make the pL sites thus 
selected unavailable for deposition for times < T < T s . A "sleeping time" 
T s for fraction p of lattice sites (grouped in dimers) in monolayer deposition 
supposedly will model effect of disallowed overhangs over gaps of size larger 
than 1 in the lower layer on the multilayer deposition in layer n provided we 
loosely identify T s oc n. Indeed, our multilayer data suggest that times needed 
to build up the nth layer coverage grow linearly with n. For instance times 
Ti/2 defined via 9 (Ti/ 2 ) = \9 n {oo), grow according to Ti/ 2 — rn, where the 
coefficient r is of order 1. After time T s all the blocked sites are released 
and can be occupied in subsequent deposition attempts. We find that the 
variation of the jamming coverage cx T~^, <j> ~ 1/3. 
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2. Continuum limit in RSA 

We consider [1] the deposition of objects of size I on a ID substrate of size L. 
R is the rate of random deposition attempts per unit time and volume. The 
lattice approximation is introduced by choosing the cubic mesh size b = l/k. 
The lattice deposition is defined by requiring that the objects of size I can only 
deposit in sites consisting of k lattice units. The late stage of the deposition 
(after time r) in continuum can be described as filling up of voids small enough 
to accommodate only one depositing object. At this time r, the density of 
those small gaps (of various sizes) will be p. For lattice models, a similar 
picture applies for k » RIt. Typical small gaps can be assumed [4] of 
sizes k + n (n = 0, 1, . . . , k — 1) with density p/k at time r, and will be 
filled up at the rate Rb(n + 1) per unit time. We will consider t » r so 
that no new small gaps are created by the filling up of large gaps. Then 
the density O of each type of the small gaps will have the time dependence 
fi(n) = (p/k) exp(—Rb(n + l)(t — r)). In each deposition event with rate 
RbQ(n)(n+l) per unit time the coverage is increased by (l/L). Thus, we have 
dO/dt ~ En=o {Rblp/k)(n + 1) exp{-i?fr(n + l)(i - r)}. After integration we 
get the asymptotic (t » r) estimate, generalized to D dimensions [1]: 



We consider some special limits. For k fixed, the "lattice" long time behavior 
sets in for Rl D t » k D . In this limit the rij = term in the sums in 
eq.(2) dominates: 6{t) w 0(oo) - (pl D /k D )exp(-Rl D t/k D ). Thus, the time 
decay constant increases as k D , The continuum limit of eq.(2) is obtained for 
k D >> Rl D t. In this limit one can convert the sums to integrals. Recall that 
all the expressions here apply only for t >> t and k D >> RI d t, where RI d t 
is a fixed quantity of order 1. Thus, the large-A; and large-t conditions are 
simply k » 1 and t » l/(Rl D ). The latter condition allows us to evaluate 
the integrals asymptotically, to the leading order for large t, which yields 



6{t) w 0(oo) - {p/(D- 1)! Rt) [In (Rl D t)] D \ The asymptotic (Int) ' 1 ^ 1 



law was derived in Ref. [4] for the continuum deposition of cubic objects. 
We evaluated 9k(oo) numerically for system sizes L/k = 200. Our data sug- 
gest a fit of the form #fc(oo) = #00(00) + (A\/k) + (A 2 /k 2 ) + By standard 

manipulations to cancel the leading l/k term, followed by a further extrap- 
olation to k — > 00, we arrived at the estimate #00(00) = 0.5620 ± 0.0002. 
The errors are small enough to rule out the conjecture of Palasti [5], and its 
generalization for finite k, which state that the jamming coverages for the 2D 
(k x k) oriented squares are equal to the squared jamming coverages of the 
corresponding ID fc-mer models. The latter are known exactly [2]. 
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3. Effects of diffusional relaxation 

In this section we report numerical studies [1] of the effects of diffusion on 
RSA in ID and 2D. In the deposition of fc-mers on a ID linear lattice, holes of 
k — 1 sites or less cannot be reached by deposition. Diffusion of the deposited 
objects can combine small holes to form larger landing sites accessible to 
further deposition attempts leading to a fully covered lattice at large times. 
For large times the holes are predominantly single-site vacancies which hop 
due to fc-mer diffusion. They must be brought together in groups of k to be 
covered by a depositing fc-mer. If the deposition rate is small, the fc-site holes 
may be broken again by diffusion before a successful deposition attempt. Thus 
the process of fc-mer deposition with diffusion will reach its asymptotic large- 
time behavior when most of the empty space is in single-site vacancies. The 
approach of the coverage to 1 for large times will then be related to the reaction 
kA — > inert with partial reaction probability on each encounter of k diffusing 
particles A. Scaling arguments indicate that the particle density for k > 3 will 
follow the mean-field law oc t _1 /( fe_1 ) for large times, with possible logarithmic 
corrections for k = 3 (borderline). This corresponds to 1 — 0(t-large) = 
t -1 /( fe-1 ) in deposition. Now we study the effect of diffusional relaxation in 
ID dimer deposition. At each Monte Carlo step a pair of adjacent sites on 
a linear lattice (L = 2000) is chosen at random. Deposition is attempted 
with probability p or diffusion otherwise with equal probability to move one 
lattice spacing to the left or right. The time step AT = 1 corresponds to L 
deposition/diffusion-attempt Monte Carlo steps. We define the time variable 
t = pT. Our Monte Carlo results were obtained for p = 0.9, 0.8, 0.5, 0.2. The 
coverage increases monotonically with (1 — p)/p at fixed pT. For p < 1 we 
obtain 0(oo) = 1, for p = 1, 9{T) = 1-exp {-2 [1 - exp(-T)]} -> l-e~ 2 < 1. 
The convergence of (1 — 9) to the limiting value at t — pT — oo is exponential 
without diffusion. Small diffusional rates lead to the asymptotically ~ i -1 / 2 
convergence to 9(p < l,t = oo) = 1. For faster diffusion, the onset of the 
limiting behavior is preceded by the region of ~ t^ 1 behavior followed by a 
crossover to <~ t -1 / 2 for larger times. In the cases k = 3,4 the large time 
results arc roughly consistent with the mean-field relation. For all p values 
studied the void area is dominated by the single vacancies precisely in the 
regime where the mean-field law sets in. For fc = 2 the single-site vacancies 
take over for t>2. For fast diffusion there follows a long crossover region from 
the initially mean-field to the asymptotically fluctuation behavior. 

Now we report on collective effects in RSA of diffusing hard squares. In 
each Monte Carlo trial of our simulation on a L x L square lattice, a site is 
chosen at random with deposition probability p. Only if the chosen site and 
its four nearest-neighbor sites are all empty the deposition is performed. A 
diffusion move by one lattice spacing is made if the targeted new site and its 
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nearest neighbors are all empty. Numerical estimates were obtained for the 
coverage and the "susceptibility" \ = L 2 [(m 2 ) — (|m|) 2 ], where the average 
( ) is over independent runs. The order parameter was defined by assigning 
"spin" values +1 to particles on one of the sublattices and —1 on another 
sublattice. The effective domain size, £(T), was defined by £ = 2L v r {ri 2 }. 
For p = 1, the approach of 9(T) to the jamming coverage 9(oo) ~ 0.728 < 1 
is exponentially fast. With diffusion, one can always reach the full coverage 
9(<x>) = 1. However, the approach to the full coverage is slow, power-law. Here 
the coverage growth mechanism for large times is due to interfacial dynamics. 
The void space at late times consists of domain walls separating spin-up and 
spin-down ordered regions. Since a typical domain has area ~ £ 2 {T) and 
boundary <~ £(T), we anticipate that for large times 1 — 9(T) oc £~ 1 (T). We 
found that the data roughly fit the power law, 1 - 9{T) oc T" 1 / 2 , for T > 10 3 . 
Thus, the RSA quantity 1 — 9(T) behaves analogously to the energy excess 
in equilibrium domain growth problems. The "susceptibility" \ f° r a given 
finite size L has a peak and then decreases to zero, indicating long-range 
order for large T. The peak location seems size-dependent, at T pea k oc L 2 . 
Since finite-size effects set in for £(T) ~ L, which given the "bulk" power law 
£{T) ~ T 1 / 2 leads precisely to the criterion T <~ L 2 , we expect this maximum 
in fluctuations to be a manifestation of the ordering process at high densities. 
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